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ABSTRACT 



We provide a simple theoretical model for the quasar luminosity function at high red- 
shifts that naturally reproduces the statistical properties of the luminous SDSS quasar 
sample at redshifts z ~ 4.3 and z > 5.7. Our model is based on the assumptions that 
quasar emission is triggered by galaxy mergers, and that the black hole mass is propor- 
tional to a power-law in the circular velocity of the host galactic halo, v c . We assume 
that quasars shine at their Eddington luminosity over a time proportional to the mass 
ratio between the small and final galaxies in the merger. This simple model fits the 
quasar luminosity function at z ~ 2-3, reproduces the normalization and logarithmic 
slope (/3 ~ —2.58) at z ~ 4.3, explains the space density of bright SDSS quasars at 
z ~ 6.0, reproduces the black hole - halo mass relation for dormant black holes in the 
local universe, and matches the estimated duty cycle of quasar activity (~ 10 years) in 
Lyman-break galaxies at z ~ 3. Based on the derived luminosity function we predict the 
resulting gravitational lensing rates for high redshift quasars. The lens fractions in the 
SDSS samples are predicted to be ~ 2% at z ~ 4.3 and ~ 10% at z > 5.7. Interestingly, 
the limiting quasar luminosity in our best-fit relation L cx v^/G, scales as the binding 
energy of the host galaxy divided by its dynamical time, implying that feedback is the 
mechanism that regulates black hole growth in galactic potential wells. 

Subject headings: Quasars: luminosity function - Gravitational lenses: lens statistics 

1. Introduction 

While the quasar luminosity function has been studied extensively at redshifts below z ~ 3 (e.g. 
Boyle, Shanks & Peterson 1988; Hartwick & Schade 1990; Pei 1995), the Sloan Digital Sky Survey 
(SDSS; Fukugita et al. 1996; Gunn et al. 1998; York et al. 2001) has in recent years substantially 
increased the number of quasars known at z > 3.5 (Fan, Strauss et al. 2001a, b; Schneider et 
al. 2001). Two samples of very high redshift SDSS quasars have been presented to date. The first 
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of these is a sample of 39 luminous quasars with redshifts in the range 3.6 < z < 5.0 and a median 
of 4.3. The second sample consists of 4 quasars at z > 5.7, including the quasar with the highest 
known redshift (z = 6.28). These samples begin to sketch out the luminosity function of quasars 
at z > 4 and are very important for studies of the ionizing background radiation and of quasar 
evolution around the epoch of reionization (Haiman & Loeb 1998). Although the logarithmic slope 
of the luminosity function for bright quasars admits a universal value of ~ —3.5 at redshifts z < 3, 
Fan et al. (2000a) find a shallower logarithmic slope of ~ —2.5 for their luminous quasars at z > 3.5. 
At first sight this is a puzzling result because the exponential tail of the Press-Schechter (1974) 
mass function describing the space density of the massive host galaxies in which bright quasars are 
thought to reside, becomes progressively steeper with increasing redshift. 

There are several published models of varying complexity to describe the observed evolution 
of the quasar luminosity function within hierarchical structure formation theory (Efstathiou & 
Rees 1988; Small & Blandford 1992; Haehnelt & Rees 1993; Haiman & Loeb 1998; Haehnelt 
Natarajan & Rees 1998; Kauffmann & Haehnelt 2000). In these models, the luminosities and 
lifetimes of the quasars are governed by the black hole mass and the supply of cold accreting gas. 
The abundance and evolution of the supermassive black holes assumed to power the quasars, are 
linked to the evolution of the mass-function of galactic halos. The density of bright quasars is 
observed to decline rapidly with redshift below z ~ 2 (Boyle, Shanks & Peterson 1988; Hartwick & 
Schade 1990). This decline has been explained (Kauffman Sz Haehnelt 2000; see also Cavaliere & 
Vittorini 2000) in terms of decreases in the merger rate and the availability of cold gas to fuel the 
black holes. In addition, the density of bright quasars is also observed to decline beyond a redshift 
of z ~ 2 (Warren, Hewett & Osmer 1994; Schmidt, Schneider & Gunn 1995; Kennefick, Djorgovski 
&; de Carvalho 1995). Under the assumption that quasar activity is related to the formation rate 
of halos, model luminosity functions have been constructed that successfully describe this decline 
and the shape of the luminosity function between z ~ 2 and z ~ 3.5 (Efstathiou & Rees 1988; 
Haehnelt & Rees 1993; Haiman Sz Loeb 1998; Haehnelt, Natarajan & Rees 1999; Haiman, Madau 
& Loeb 1999). These models predict the luminosity function at still higher redshifts. With the 
recent influx of results on the quasar luminosity at redshifts above z ~ 4, it is interesting to revisit 
the question of modeling the evolution of the quasar luminosity function at high redshifts. In this 
paper, we demonstrate that a very simple model can explain the new results on the high redshift 
quasar luminosity function. 

Our model for the evolution of the quasar luminosity function is based on the Press-Schechter (1974) 
mass function and halo merger rates computed with the excursion set formalism (Bond et al. 1991; 
Lacy & Cole 1993). We show that this model can reproduce all of the known properties of the 
luminosity function above z ~ 2. We begin in §2 by describing an approximate version of our model 
(Haiman & Loeb 1998; hereafter HL98), followed by the full discussion based on the halo merger 
rates. In §3 and §4 we discuss the data on the quasar luminosity function and its comparison with 
our model. Finally we discuss the implications of our model for gravitational lensing in high-redshift 
samples of quasars in §5, and present our main conclusions in §6. Throughout the paper we assume 
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density parameters of £l m = 0.35 in matter, f^A = 0.65 in a cosmological constant, and a Hubble 
constant of Hq = 65 kms _1 Mpc _1 . For calculations of the Press-Schechter (1974) mass function 
we assume a primordial power-spectrum with n = 1 and the fitting formula to the exact transfer 
function of Bardeen et al. (1986). 



2. Theoretical Luminosity Functions 

We compare two models for the quasar luminosity function that are based on the Press- 
Schechter (1974) mass function of dark matter halos, and the halo merger rates from the excursion 
set formalism of Bond et al. (1991) and Lacy & Cole (1993). We begin with the approximate model 
of HL98 which is based on the halo formation rate. We then describe our more complete model 
based on halo merger rates. We will demonstrate in §4.1 that only this latter model reproduces all 
the known features of the quasar luminosity function above z ~ 2, as well as the local black hole - 
host halo mass relation and new constraints on the quasar duty cycle. 



2.1. An Approximate Model for the Quasar Luminosity Function Based on the 

Halo Formation Rate 

HL98 discussed an approximate model for the high redshift quasar luminosity function by 
supposing that quasars are associated with newly formed halos. Their model assumes that each 
galactic halo hosts a black hole with a mass (Mbh) which is a constant fraction (e) of the host halo 
mass (Mh a i ), and that the black hole shines at its Eddington rate with a universal light-curve, 
f(t). HL98 used the derivative of the Press-Schechter mass function to approximate the formation 
rate of halos, arguing that mergers are rare at high redshift so that the negative contribution to 
the derivative is negligible. In this model the quasar luminosity as a function of time is 

L B (t) = M hh f(t) = eM hSLlo f(t) for M halo > M min , (1) 

where M m - m ~ 10 8 Mq[(1 + z)/10] _ 2 is the minimum halo mass inside which a black hole can form. 
This lower mass limit corresponds to the virial temperature below which atomic cooling is not 
effective in allowing the gas to sink to the center (Barkana & Loeb 2001). Following HL98 we write 
the resulting quasar luminosity function, defined as the comoving number density of quasars having 
rest frame B-band luminosities between Lb and Lb + ALb and redshifts between z and z + Az, as 

*hl(L b , z) = j°° dz' dM hh d Ub] h 6[L B - M hh f(t z - t% (2) 

Jz JeM min dM hh dZ 

where t z and t' are the time at redshifts z and z', and d ^" h b ^ z ; is the change in the co-moving 
number density of black holes between z' and z' + dz' . Integrating over we get 
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where a = (1 + z)~ l is the scale factor at a redshift z. HL98 used an exponential with an e-folding 
time idc,o f° r the universal light curve /. We find that the equivalent result for ^hl(-^B) z) is 
obtained more simply by assuming a step function for /, namely 

f(t) = ^fe(t-t dc , ), (4) 

where -£/Edd,B = 5.7 x 10 3 (M^/Mq) is the Eddington luminosity in B-band solar luminosity units 
for the median quasar spectrum (Elvis et al. 1994), and Q(t) is the Heaviside step function. We 
then find 

1 a= T+z dz M bh (i 2 nbh 



*hl(£b,*) = I 



da 

i daL EddtB dM hh dz 



(5) 



M bh = 



f(tz-t) 



1 ' z forn 

given the formation redshift for the halo, Z{ OTm . Assuming that tdc,o ^ H^ 1 (z), (a — a m - m ) « 
(da/dt)tdc,o, and relating the space density of black holes to the Press-Schechter (1974) mass func- 
tion through 



d 2 n hh 1 d 2 n ps 

(6) 



dM^hdz e dM\y & \ dz ' 
we obtain the following simple expression for the quasar luminosity function: 



\[l (T \ ^dc,0 &m , 1 1 d 2 Ups 

HLl B,Z) ~ jy- 1 V a 5 a 2 5.7 x 10 3 e dM h&lo dz 
We will compare this model to observations in §4.2. 
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2.2. A Model for the Quasar Luminosity Function Based on the Halo Merger Rate 

The HL98 model associates quasar activity with halo formation. The assumption is that the 
halo merger rate is rarer than the halo formation rate at high redshift. However, at low redshifts 
the formation rate becomes smaller than the merger rate, and the derivative of the Press-Schechter 
mass function on which the model is based, becomes negative. Quasar activity at low redshift can 
be explained in relation to halo merger activity (Carlberg 1990; Kauffmann & Haehnelt 2000). In 
this section we discuss a simple model that associates quasar activity at high redshift with the 
merger rate of halos. 

The merger rate of halos was computed by Lacy & Cole (1993) based on the excursion 
set formalism of Bond et al. (1991). We compute the number of halos having masses between 
AMh a io and AMh a i + dAMhaio that accrete onto a halo of mass Mh a io — AJV4 a i per unit time 
. Thus, the number of merger events involving the accretion of a halo of mass 
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given by the product 
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AMhaio by a halo of mass Mh a i — AMhaio per unit time per comoving volume at a redshift z is 
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M=M halo -AM halo dAM halo dt 
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We assume that the mass of the central black hole scales as a power-law with the circular 
velocity (v c ) of the halo, 

M bh oc v2. (9) 



The circular velocity of a halo of mass Mh a i at redshift z can be written as (Barkana & Loeb 2001), 
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where /i = (#o/100 km s 1 Mpc 1 ), A c = 187r 2 + 82c? — 39d 2 is the final overdensity relative to the 
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and fit for two free parameters e Q and 7 that describe the evolution and slope of the black hole 
- halo mass relation. We assume that the black holes coalesce upon halo merger (Kauffmann & 
Haehnelt 2000) and find the number of black holes of mass between AM D h and AM D h + dAM D h 
that merge with black holes of mass Mbh — AMbh per unit time: 
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(12) 
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A linear relation (7 = 3) between black hole mass and halo mass is a natural consequence of black 
hole growth that is dominated by coalescence with no gas accretion (Haehnelt et al. 1998). Values 
of 7 > 1 result from significant gas accretion during the active quasar phase as deduced by Yu 
&; Tremaine (2002). We assume that after a merger, a fraction of the cold gas from the accreted 
halo is driven onto the central black hole of mass Mbh (Mihos k, Hernquist 1994; Hernquist & 
Mihos 1995). If the quasar shines at the Eddington rate of the black hole in the merger product 
(see, e.g. Yu & Tremaine 2002), then the cold gas from the small accreted halo (which makes up 
the new fuel reservoir) will run out in a time approximately proportional to 
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(13) 



where AMb ary0 n is the mass in baryons within the accreted halo 2 . Thus we postulate that after a 
merger the quasar shines at the Eddington rate corresponding to the merger product, for a time 
proportional to both the ratio between the masses of the accreted and initial halo and the baryon 
fraction of the accreted halo. 

In order to examine whether the baryon mass fraction in the accreted halo would be affected by 
the photo-ionization heating of the intergalactic medium (IGM) after reionization, we have solved 



2 Equation (13) can be derived as a Taylor expansion for minor mergers. For simplicity we apply it to all mergers. 
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the linear growth factors D^ in and D b for the dark matter and baryons on a spatial comoving scale 
R = (3AMh a io/47r / o m ) 1 ^ 3 , where p m is the average matter density of the universe today. These 
growth factors obey the following set of coupled differential equations for the dark matter and 
baryon overdensities 5dm and 5 b in the linear regime (e.g. Barkana & Loeb 2001): 



<5dm + 2F<Wi = -^H 2 [£l b (z)5 b + fidmO)<W] 
5b + 2H5 h = Ih 2 {[n b (z)5 b + fW^dm]} - 



fj,m p 



- (-) p) 



^ + 3/?(5 b - <y b) i) 



(14) 



Here T is the gas temperature, fi the mean molecular weight (~ 0.69 for ionized IGM), fee is 
Boltzmans constant, and k is the co-moving wave number for the mode of interest. The subscripts 
i refer to quantities at an initial reference time and (3 = or 1 for adiabatic or isothermal evolution. 
We assume reionization at z TC \ OIi = 7. Before reionization we assume 7] = 0; the baryon and dark- 
matter over-densities undergo the same evolution during this period. Following reionization, the 
IGM is assumed to be heated to 10 K and to evolve isothermally (/? = 1). The evolution of D b 
is described by the superposition of solutions to equation (14) weighted by the Fourier transform 
of the spatial tophat window function of width R (Medvigy & Loeb 2001). Figure 1 shows the 
evolution of D b with redshift. The curves shown correspond (from bottom to top) to halo masses 
of M ha i = 10 6 , 10 7 , 10 8 , 10 9 , 10 10 , 10 11 , 10 12 , 10 13 and 1O 14 M . The dark matter growth factor 
corresponds to the upper envelope in this figure. Prior to reionization D b follows -Ddm f° r an halo 
masses. Reinonizaton heats the IGM and eliminates gas accretion onto small halos. The baryon 
overdensity in halos larger than 10 10 M Q is unaffected by reionization. Thus, we find that the effect 
of reionization on the luminosity function is only apparent for quasars fainter than ~ 1O 11 L , 
which is well below the detection threshold of existing quasar surveys at z > 1, such as SDSS. 
We therefore assume a fixed baryon mass fraction in objects relevant for the luminosity function of 
bright quasars at high redshifts. 

The quasar light curve and luminosity function may therefore be written as 
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in analogy with equations (4) and (5), respectively. Integrating over Mbh we find 
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Assuming td c ,o ^ H 1 (z), we then obtain a — a mm = ^td c ,o Z ^ fhal ° , yielding the luminosity function 
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where M ha i = L B /5.7 x 10 3 e. 
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(18) 



The value of tdc,o can be related to the duty cycle of quasars i^c through consideration of the 
number of mergers during a Hubble time H~ l (z), 
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We will compare the luminosity function predicted by this model with observations in §4.1. 



3. The Observed Luminosity Function 

The standard double power-law luminosity function (Boyle, Shanks & Peterson 1988; Pei 1995) 

= wwi-y (20> 

provides a good representation of the observed quasar luminosity function at redshifts z < 3. At 
the faint end of the luminosity function, the slope is 0\ = 1.51, while at the bright end j\ = 3.43. 
Moreover, all dependence on redshift (at z < 3) is in the break luminosity indicating pure 
luminosity evolution. At low redshift the break luminosity evolves as a powerdaw in redshift, and 
the space density of bright quasars increases with redshift. However, surveys at higher redshift 
show a decline in the space density of bright quasars beyond z ~ 3 (Warren, Hewett & Osmer 1994; 
Schmidt, Schneider & Gunn 1995; Kennefick, Djorgovski & de Carvalho 1995). 

While quasar evolution below z ~ 3 is described by pure luminosity evolution, Fan et al. (2001a) 
found from their sample of SDSS quasars at z ~ 4.3 that the slope at the bright end of the luminosity 
function has evolved from the z < 3 value of (3h ~ —3.5 to (5^ ~ —2.5. This result is supported 
by the analysis of Schmidt, Schneider & Gunn (1995). Fan et al. (2001a) also find an evolution of 
space density with redshift for bright quasars of ^ oc l0 -a52: between z ~ 3.5 and z ~ 5. The space 
density of bright quasars measured by Fan et al. (2001b) at z ~ 6 agrees with the extrapolation of 
this evolution. 

We seek to model the observed features of the luminosity function at redshifts higher than the 
peak in quasar evolution at z ~ 2. In the following section we will discuss both our model (§2.2), 
and the HL98 model in light of the recently acquired data on quasars at high redshifts. 
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4. Comparison of Observed and Model Luminosity Functions 

In this section we compute model luminosity functions at different redshifts and compare them 
with the observed quasar luminosity function. Figure 2 shows the observed luminosity function at 
z = 0.1, 1.0, 2.3, 3.0, 4.0, z ~ 4.3 and z ~ 6.0. Below z ~ 4.3, the data shown is from the summary 
by Pei (1995; based on the compilations of Hartwick & Schade 1990). We also show the best fit 
empirical luminosity function at z ~ 4.3 from Fan et al. (2001a) with a vertical bar to denote the 
quoted error in the normalization. At z ~ 6 the space density was inferred from Fan et al. (2001b), 
and the error bar includes uncertainty in the spectral index. Comparisons of this data with the 
halo merger rate model (§4.1) and the HL98 model (§4.2) are presented below. 

4.1. The Merger Rate Model 

Luminosity functions calculated from the halo merger rate model described in §2.2 are plotted 
as the solid lines in Figure 2 at redshifts z = 0.1, 1.0, 2.3, 3.0, 4.0, z ~ 4.3 and z ~ 6.0. In addition, 
the lower right panel shows predicted luminosity functions at z = 8.0 and z = 10.0. Values of 
e G = 10~ 5 ' 2 , 7 = 5 and tdc,o = 10 6 ' 3 yr yield a good fit to the data at all z > 2. The parameters 
in the black hole - halo mass relation bear a striking resemblance to those in the relation found 
by Ferrarese (2002) for the local universe, which in our notation correspond to e G = 10 -51 and 
7 = 4.71. 

The value of tdc,o is related to the duty cycle for different quasar luminosities and redshifts 
through equation (19). The duty cycle is plotted in Figure 5. For a given halo mass, t& c is longer 
at higher redshift. Steidel et al. (2002) estimated the lifetime of bright quasar activity in a large 
statistical sample of Lyman-break galaxies at z ~ 3 to be ~ 10 7 yr, which is surprisingly close to 
our expected value of t& c = 10 6 ' 9 yr. Above z ~ 3 we find tdc ~ 10 7 yr, comparable to the value 
inferred by comparing the local black hole density with the quasar luminosity function (e.g. Yu & 
Tremaine 2002). Note that while tdc becomes a substantial fraction of the age of the universe at 
the highest redshifts considered, individual episodes of quasar activity are shorter than idc,o> which 
is an order of magnitude smaller. 

Two results from the comparison between model and data stand out. First, the halo merger 
model reproduces the shallow slope for the z ~ 4.3 quasars observed by Fan et al. (2001a), and 
the space density of quasars at z ~ 6.0 measured by Fan et al. (2001b). A detailed comparison 
between the data and the model at z ~ 4.3 is shown in Figure 3. The left hand panel shows the 
luminosity function in the region measured by Fan et al. (2001a). The straight line with the error 
bar represents the empirical luminosity function (best fit slope of f3 = —2.58) and the uncertainty 
in its normalization. The upper solid line is the luminosity function of the merger rate model. 
The theoretical luminosity function slope is consistent with that of the observed quasars over the 
measured luminosity range. This is shown more clearly on the right hand panel of Figure 3 where 
the derivative of the model luminosity function is shown in comparison with the empirical best 
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fit slope (solid line) and the slope ±1(7/3 away (dashed lines). The model is consistent within the 
quoted uncertainty over the measured luminosity range, particularly at the fainter end. Our model 
luminosity function is consistent with the observed lower limit on the slope at z = 6.0, which Fan 
et al. (2000b) quote as /?h < 3.9 (95%). However the gravitational lensing rate may be very high in 
this sample (Wyithe & Loeb 2002a,b) . While the average slope of the luminosity function cannot 
be changed significantly by lensing, we might easily be misled by the magnification of 1 or 2 quasars 
in the limited current sample of 4. Furthermore, lensing might increase the apparent space density 
of quasars by a factor of up to 2 (Wyithe & Loeb 2002b). The predicted luminosity functions show 
that the slope of the bright end becomes significantly steeper at redshifts higher than z ~ 6 as 
expected from the shape of the Press-Schechter (1974) mass function. 

The fit to the luminosity function discussed above has 7 = 5. An equally good fit can be 
obtained for 7 = 4 and e Q = 10 -4 ' 5 . However this value of e Q is nearly an order of magnitude 
larger than observed (Ferrarese 2002), and the implied duty cycle is nearly an order of magnitude 
too small (Yu & Tremaine 2002; Steidel et al. 2002). Hence, the complete set of observational 
constraints favors a value of 7 = 5. Ferrarese & Merritt (2002) find a relation between Mbh 
and bulge velocity dispersion of Mbh °^ u4.72io.36 while Tremaine et al. (2002, and references 
therein) find Mbh °^ a"f' 02=l=0 ' 32 . Tremaine et al. (2002) attribute this disagreement primarily to 
the differences in the velocity dispersion measurements used by the two groups. Unfortunately, 
the highly uncertain nonlinear scaling of v c oc a ^ 88±0 - 17 (Ferrarese 2002) allows the results of both 
groups to be consistent with our favored value of 7 = 5 (which in turn translates to a mean slope 
of Mbh « 0"b' 4 ' in between the slopes of the two groups). 

The upper left hand panel of Figure 2 indicates that the model luminosity function provides 
a poorer fit to the data at redshifts z < 2. The density of relatively faint [Lb < 1O 13 L at z ~ 1.0 
and Lb < 1O 12 L at z ~ 0.1 (dashed lines)] quasars is reasonably reproduced, particularly at the 
lowest redshift considered. The model also predicts a decline in the number of bright quasars at low 
redshifts. On the other hand, the slope of the model luminosity function is too shallow at the higher 
luminosities to explain the data. A simple observational explanation for this phenomenon might 
be the large drop in the average cosmic density of cold galactic gas from z ~ 2-3 to the present 
day (Storrie-Lombardi, McMahon & Irwin 1996). Furthermore, the inability of the gas to cool 
following the epoch of group and X-ray cluster formation would also contribute to the rapid decline 
in the luminosity function at low redshifts (Cavaliere & Vittorini 2000). Thus, while the merger 
rate predicts a decrease in the space density of quasars at low redshifts, the more complicated low 
redshift environments require additional modeling (e.g. Kauffmann & Haehnelt 2000) to reproduce 
the observed luminosity function shape. 
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4.2. The Halo Formation Model 

We find that the HL98 model (§2.1) requires 3 best-fit values of e = 1CT 3 - 1 and tdc,o = 10 5 ' yr. 
The resulting luminosity function is shown as the dotted lines in Figures 2 and 3. The HL98 model 
and the halo merger model produce similar luminosity functions at z > 2. However, the predicted 
duty cycle is too low by two orders of magnitude, and the value of e produced by this fit is too high 
by 2 orders of magnitude to be consistent with observations. 

Figures 3 and 4 indicate that the HL98 model produces a somewhat less consistent normaliza- 
tion at z ~ 4.3 and a steeper slope at z ~ 6.0 when compared with the halo merger model. At low 
redshifts the merger rate is higher than the collapse rate for halos in the mass range of interest. In 
fact, the derivative of the Press-Schechter (1974) mass function becomes negative and the model 
predicts a luminosity function with a sharp lower cutoff in luminosity (in the implementation of 
HL98 the smooth function, f(t), results in a softer cutoff). This is seen in the upper left panel of 
Figure 2. If quasar activity was associated solely with the halo formation rate, then only quasars 
residing within the very largest halos would be observed at low redshift. 

5. Rates of Multiple-Image Lensing 

The fraction of quasars that are gravitationally lensed by foreground galaxies is very sensitive 
to the underlying luminosity function due to the potentially large magnification bias (Turner 1980; 
Turner, Ostriker & Gott 1984). At low redshifts (z < 3) the bright end of the luminosity function 
has a universal slope of —3.43. However at z ~ 4.3 Fan et al. (2001a) observed a slope of —2.58. 
Wyithe & Loeb (2002a,b) showed that empirical extrapolations of the low redshift quasar luminosity 
function to high redshifts yielded very different results for the gravitational lensing rate if one were 
to adopt a bright end slope that is —3.43 at all redshifts or a slope that is —2.58 above z ~ 3. 
Since the halo merger model described in §2.2 reproduces all measured properties of the quasar 
luminosity function at high redshifts, we use it to compute the magnification bias for multiple- 
image gravitational lenses. The magnification bias can then be used to find the fraction of all 
quasars (at different redshifts) that are multiply imaged. The lensing fraction for the SDSS high 
redshift samples was discussed at length in Wyithe & Loeb (2002a,b). We therefore give only a 
short summary here, and present results for the lens fraction obtained with the new luminosity 
function. 

The probability distribution of point source magnification by a singular, spherical isothermal 
lens in terms of the sum of the magnifications of multiple images, fi, is {dP/dfi) = 8/fi 3 for n > 2. 
The bias factor for multiple imaging in a flux limited sample brighter than (absolute luminosity) 



3 HL98 quote tdc,o = 10 5 ' 8 yr, larger than our value of to = 10 5 . Haehnelt, Natarajan & Rees (1998) also found a 
lower value tdcfi in their model with constant e. 



- 11 - 



Lb lim is therefore 

jf^c,.../,^*(^> .... 
() ~ JZ w.>) ' ( ' 

This magnification bias is plotted as a function of the limiting (absolute) luminosity in Figure 6. 
Biases at redshifts of z = 2.3, z = 4.3 and z = 6.0 are shown, as well as predictions for z ~ 8.0 and 
z ~ 10.0 (bottom to top). The bias becomes very large at high redshifts and for bright limiting 
magnitudes. 

The resulting fraction of multiple images expected in the flux limited sample is 

FiLa^z)* BT ™\ (22) 

In this expression r mu it is the multiple image optical depth, i.e. the source plane probability 
that a source will be multiply imaged for the observer. From Wyithe & Loeb (2002b) we find 
r muit = 0.0040 for the SDSS quasars at z ~ 4.3 and r mu i t = 0.0059 for the SDSS quasars at z ~ 6. 
In addition, we find r mu ] t = 0.0020, r mu i t = 0.0076 and r mu i t = 0.0091 (considering all galaxies as 
potential lenses) at z = 2.3, z = 8.0 and z = 10.0, respectively. The lens fraction F(L-Bfi m , z) is 
plotted on the right hand panel of Figure 6 for z = 2.3, 4.3, 6.0, 8.0 and 10.0 (bottom to top). 
While at low redshifts the lensing fraction is less than one percent, at higher redshifts the fraction 
could be much higher. The vertical dashed lines correspond to the limiting absolute luminosities 
of the SDSS quasar samples at z ~ 4.3 (left) and z ~ 6 (right). These limits suggest lens fractions 
of ~ 2% in the SDSS sample at z ~ 4.3 and ~ 10% in the SDSS sample at z > 5.7. Note that 
the fractions at z ~ 4.3 and z ~ 6.0 are consistent with the value obtained for the flat luminosity 
function assumed in Wyithe & Loeb (2002b). At even higher redshifts the lens fraction approaches 
unity. The luminosity function at z ~ 10 predicts that quasars having luminosities brighter than a 
limit that corresponds to the same co-moving space density as the SDSS quasars at z ~ 6.0 have a 
lens fraction of F ~ 0.5. 



6. Summary and Discussion 

We predict the quasar luminosity function under the hypothesis that quasar activity is asso- 
ciated with galaxy mergers. Galactic halos are assumed to each host a black hole having a mass 
determined from a power-law relation with halo circular velocity, Mbh oc v2- Upon merger we 
assume that the gas from the accreted halo is driven to the center, and that the quasar then shines 
at its Eddington rate for a short period of time, proportional to the (smaller than 1/2) mass ratio 
between the small and final halos in the merger. Our model has 3 free parameters but at least the 
following 5 constraints: the luminosity function slope and normalization at z ~ 2.5, the evolution 
of the luminosity function slope and normalization, and the curvature of the luminosity function. 
Our model is simplified by the fact that the effect of reionization on the luminosity function is 
limited to luminosities much fainter than those detectable in current quasar surveys. 
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This very simple model reproduces all known properties of the luminosity function above z ~ 2, 
including the recently measured luminosity function slope at z ~ 4.3 (Fan et al. 2001a) and the 
space density of quasars at z ~ 6.0 (Fan et al. 2001b). The model suggests a black hole - halo 
mass relation, Mbh oc consistent with that found by Ferrarese (2002) at z = 0, and a duty cycle 
of ~ 10 7 yr consistent with previous determinations (Yu & Tremaine 2002; Steidel et al. 2002). 
Using the derived luminosity function we predict the resulting lensing rates for the highest redshift 
samples. The recently published flux limited samples of SDSS quasars at z ~ 4.3 and z > 5.7 are 
predicted to have multiple image fractions of ~ 2% and ~ 10%, respectively. 

The model luminosity function computed using the halo merger rates is consistent with the 
density of faint quasars at low redshifts, and predicts a decline in the density of bright quasars 
towards low redshifts, but predicts a luminosity function slope that is shallower than observed. As 
has been noted previously, a more careful analysis (e.g. Kauffmann & Haehnelt 2000) is required 
to understand the evolution of the quasar luminosity function at low redshifts. Nevertheless, the 
consistency of the faint quasar normalization stands in contrast to models that associate quasar 
activity with the halo formation rate, as those predict no faint quasars at low redshifts. Thus one 
is led to conclude that mergers are the likely trigger for quasar activity at all redshifts. 

The Mbh-^c relation that provides the best fit to the data may also have a simple physical 
origin. First, consider dimensional analysis. If we make the minimal assumption that the limiting 
luminosity of a quasar is only a function of the halo circular velocity v c , then the only dimensional 
parameters in the problem are v c and G (Newton's constant). The only combination of these 
parameters that has dimensions of luminosity is v^/G and so we conclude 

L E dd oc ^. (23) 

What is the physics behind this relation? The value of v^/G amounts to depositing the entire 
binding energy of a self-gravitating system during its dynamical time (so that it does not have time 
to adjust). The binding energy of a self-gravitating mass M is ~ Mv\ and the dynamical time is 
~ r/v c . The ratio between these quantities is ~ Mv\jr. Using the virial relation GM/r ~ v%, we 
get the energy deposition rate that would unbind a self-gravitating system on its dynamical time, 
~ v^/G. A quasar may therefore unbind the gas in the galaxy around it 4 if its power output is 
too large 5 (Silk & Rees 1998; Ciotti k Ostriker 2001). As the mass of the black hole increases 
and the quasar's Eddington luminosity approaches this limit, the feedback will generate a powerful 
galactic wind and terminate the accretion that feeds the quasar. Hence, the quasar phenomenon 



4 We implicitly assume that the central region of the galaxy (such as the proto-bulge) is characterized by a circular 
velocity similar to that of the galactic halo, even though it is dominated by baryons. This assumption follows naturally 
from the nearly fiat rotation curves of nearby galaxies. 

Note that the power output from the quasar may include both radiation and mechanical energy. Substantial 
outflows are inferred to exist in broad absorption line quasars and radio galaxies [see recent summary in Furlanetto 
& Loeb (2002), and references therein]. 
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may be self terminating (similarly to the formation of a proto-star). Of course, if some of the 
quasar energy escapes (due to incomplete absorption, partial covering factor or efficient cooling by 
the surrounding gas), then the feedback- limited quasar luminosity would be higher than v\jG. As 
long as the feedback from all quasars encounters self-similar conditions in different halos, we may 
write 

5 

£ E dd = /3^, (24) 

where (3 e is a constant that is related to the inferred value of e Q . Equations (10), (11), (24) and the 
best-fit value of e Q ~ 1CT 5 - 2 together with the relation L E dd = 1-4 x lO 38 (M b h/M ) erg s" 1 , yield 
(3 e ss 50. We have shown that the assumption of a constant (3 e in equation (24) leads to a luminosity 
function that describes the data over a wide range of redshifts. The constancy of (3 e (which is also 
expected from dimensional analysis) implies that feedback may indeed be the mechanism that 
regulates the growth of supermassive black holes in galactic potential wells. The similarity between 
our deduced values of 7 = 5 and e Q = 10~ 5 ' 2 at high redshifts and the values inferred for remnant 
black holes in the local universe (Ferrarese 2002; Tremaine et al. 2002), provides a strong testimony 
to this effect. 

This work was supported in part by NSF grants AST-9900877, AST-0071019 for AL. JSBW 
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operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
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z 



Fig. 1. — The growth factor (normalized to the present day -D^o) for baryons associated with 
dark matter halos of mass M h£do = 10 6 , 10 7 , 10 8 , 10 9 , 10 10 , 10 11 , 10 12 , 10 13 and 1O 14 M (bottom 
to top). The down turn of the growth factor for small halos corresponds to the reionization epoch 

at Zreion ~ 7-0. 
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Fig. 2. — Comparison of the observed and model rest frame B-band luminosity functions. The 
data at z < 4 is summarized in Pei (1995). At z ~ 4.3 and z ~ 6.0 the data is from Fan et 
al. (2001a;2001b). In the z ~ 4.3 panel, the diagonal line shows the best fit slope of —2.58 measured 
by Fan et al. (2001a), and the vertical bar shows the quoted uncertainty in the normalization. In 
the lowest redshift panel, the dashed line and the stars represent the merger model luminosity 
function and the observed luminosity function data at z ~ 0.1 respectively. The dotted lines show 
the HL98 model (e = 10 -31 , t dc ,o = 10 5 '°yr) and the solid lines are the merger model described in 
this paper (e = 10" 5 - 2 , 7 = 5, t dcfi = 10 6 - 3 yr). 
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Fig. 3. — Detailed comparison of the observed rest frame B-band quasar luminosity function with 
the models at z ~ 4.3. The left panel shows the luminosity function in the region measured by Fan 
et al. (2001a). The straight line with the error bar shows the empirical luminosity function (best fit 
slope of f3 = —2.58) and the quoted uncertainty in the normalization. The upper solid and dotted 
lines show the luminosity functions for the halo merger model (eo = 10 -5 ' 2 , 7 = 5, tdc,o = 10 6 ' 3 yr) 
and for the HL98 model (e = 10~ 31 , tdc,o = 10 5 yr), respectively. On the right panel we show the 
derivatives of the model luminosity functions. The solid and dotted lines show our model and the 
HL98 model , respectively. The empirical best fit slope (solid) and the slopes ±lcrg = ±0.23 away 
(dashed) are shown for comparison. 
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Fig. 4. — Detailed comparison of the observed rest frame B-band quasar luminosity function with 
the models at z ~ 6.0. The left panel shows the space density of quasars measured by Fan et 
al. (2001b). The solid and dotted lines show the luminosity functions for the halo merger model 
(e = 10~ 5 - 2 , 7 = 5, t dcfi = 10 6 - 3 yr) and for the HL98 model (e = 10~ 31 , t dcfi = 10 5 -°yr), 
respectively. On the right panel we show the derivatives of the model luminosity functions. The 
empirical lower limit (95%) is shown for comparison. 
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Fig. 5. — The duty cycle time t^c of quasars with luminosity Lb calculated from equation (19) at 
redshifts z = 0.1, 1.0, 2.3, 3.0, 4.0, 4.3, 6.0, 8.0 and 10.0 (bottom to top, the thick line denotes the 
value at z = 3.0 which should be compared with ~ 10 7 yr found by Steidel et al. (2002)). The duty 
cycle was calculated using our best fit value of tdc,o = 10 6 ' 3 yr per merger. 
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Fig. 6. — Left: The magnification bias as a function of absolute luminosity limit for flux limited 
samples of quasars at z = 2.3, 4.3, 6.0, 8.0 and 10.0 (bottom to top). Right: The multiple-image 
lensing rate for flux limited samples of quasars at z = 2.3, 4.3, 6.0, 8.0 and 10.0 (bottom to top). 
The left and right vertical dashed lines correspond to the absolute flux limit of the SDSS samples 
at z ~ 4.3 and z ~ 6.0, respectively. 



